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ABSTRACT 



Aim of this work is a broad-band study with INTEGRAL, Swift and XMM- 
Newton satellites of a sample of 9 blazars (7 FSRQ and 2 BL Lac) with redshift up to 
about 4. 

The spectral analysis has shown clear evidence of a flattening of the continuum to- 
wards the low energies (E < 3 keV observer frame) . This behaviour is well reproduced 
both with an absorbed power-law model (Njj ~ 10 20 -10 23 cm~ 2 in the rest-frame of 
the sources) or a broken power-law continuum model (with an energy break below 3 
keV in the observer-frame). No Compton reflection features, Fe Ka line and hump at 
high energies, have been detected, with the exception of the source IGR J22517+2218 
that shows the presence of a weak iron line. 

In this work we also investigate a possible correlation between the absorption 
column density Njj and the red-shift. 

We confirm the existence of a A#-z trend, with the higher absorption at z>2 
for a larger sample compared to previous results. The distribution of the Nh and 
the photon index T is also presented. The hard X-ray data allow us to detect highly 
absorbed sources (with Nh ^ 10 23 cm~ 2 in rest-frame of the source) characterized 



by photon index distribution peaked at harder values (F 
obtained with XMM data only (r ~ 2). 



1.4) with respect to that 
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1 INTRODUCTION 

Blazars are one of the most intriguing class of objects among 
the Active Galactic Nuclei (AGN) and include BL Lac- 
ertae objects (BL Lac) and Flat-Spectrum Radio Quasars 
(FSRQ). The intrinsic differences between these two types 
of sources are in their optical spectra and in the intensity 
of their emission: a BL Lac source does not show any op- 
tical emission lines (EW < 5A) and it is basically a low- 
power blazar (Lsoi ~ 10 46 — 10 47 erg s _1 ), whereas a FSRQ 
source shows significant emission line equivalent widths and 
corresponds to a h i gh-power blazar (LboI ~ 10 48 erg s _1 ) 
JStickel et all Il99l|: iPadovanil ll997l L In the AGN unified 
scheme ( Antonucci 1 19931 ). blazars are interpreted as radio- 
sources with a relativist ic jet aligned along the line of sight 
(|Urrv fe Padova ni 1995); in other words, the class of blazars 
is represented by Radio Loud (RL) AGN observed very close 
to the direction of the relativistic jets (< 10°). 

The Spectral Energy Distribution (SED) of blazars 
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is usually modelled with two large humps produced via 
Synchrotron (peaking in the infra-red to soft X-ray en- 
ergy band) and Inverse Compton (IC) emission, dom- 
inating the hard X-ray to gamma-ray regimes, respec- 
tively. The lower energy emission component is Synchrotron 
radiation from the jet whereas the high-energy compo- 
nent arises through the inverse Compton of soft pho- 
tons by highly relativistic electrons in the jet plasma. 
These soft photons originate from the local Synchrotron 
radiation within the jet or the nuclear optical/UV 
emission, namely Synchrotron Self-Compton (SSC) and 
External Compton (EC) compon e nts, respectively (see 
i Maraschi. Ghisellini fe Celottil Il992l: iDermer fe Schlickeiserl 
ll993l ; ISikora. Begelman fe ReesHl99€) . 

The X-ray spectrum of an AGN is usually well described 
by a power-law with a specific flux (i.e. per unit energy inter- 
val) of the form N(E) oc E~ T , where E is the energy, N(E) 
is the number of photons in units of s _1 cm -2 keV -1 and Y 
is the photon index. However, X-ray spectra of blazars show 
some deviations from the simple power-law. These spectral 
signatures appear in the soft X-ray band with a curvature (a 
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flattening or a steepening of the spectrum); in addition, the 
presence of Compton reflection components in RL objects is 
still debated. 

A clear physical interpretation of these features has not 
yet been found. In the present work, we investigate these still 
open questions in order to achieve a better understanding of 
the general properties of the blazars spectra. The state-of- 
art of the study of these important issues can be summarize 
as follows: 



(i) Two physical interpretations to the flattening of the 
primary intrinsic co ntinuum wi t h respect to a simple power- 
law exist so fa r (|Page et all 120051 ; iGalbiati et al] 120051 ; 
I Yuan et al-lbood): absorption in excess of the Galactic com- 
ponent (e.g. lCappi et all 19971: lYuan et al . 2006) and a break 
in the intrinsic contin uum (|Sikora. Begelman fe Reesl 1 19941 : 
iTavecchio et al.ll2007n . For the absorption excess, a depen- 
dence between the hydrogen column density of the absorber 
{Nh) and the rcdshi ft z has be en proposed. In fact in a for- 
mer work bv lYuan et al] (|2006l ) based on XMM observations 
of a sample of 32 RL sources, a Nh-z correlation has been 
found. 

(ii) Contrary to Radio Quiet (RQ) sources the pres- 
ence of reflection features - most notably the Fe Ka 
line and the associated Compton reflection "hump" at 
about 20 - 30 keV - from the cores of the RL AGN 
is not well established. In particular, the reprocessed 
feature s are genera lly intense and always detected in RQ 
AGN (|Mattl |2003 ). whereas in RL AGN the iron line 
and the reflection component can be absent or weak. 
X-ray observatio ns of RL AGN with ASCA, RXTE and 
Beppo SAX (e.g. Sambruna. Eracleous fe Mushotzkvl 19991: 



Eraclcous, Sambruna & Mushotzkv 2000; Grandi et al 
2001] ; iBallantvne. Ross fe Fabianl |2002| ; iGrandi fe Palumbo 



2004 ; IGrandi Malaguti fe Fiocchil 12009 ) have shown that 



some of these objects (defined as Broad Line Radio 
Galaxies) seem to have weak hard X-ray reflection features. 

In this paper we present a broad-band (0.2 — 100 keV) 
study of a sample of 9 blazars up to redshift ~ 4 observed 
with INTEGRAL, Swift and XMM-Newton satellites. The 
wide energy range covered by the instruments is well suited 
to better address their spectral behaviour. 

In Sections 2.1-2.4 we present the data set and analysis, 
while in Section 3.1 we show the details of the spectral anal- 
ysis. The results are discussed in Section 3.2 and in Section 
3.3 a study of the distribution of different spectral param- 
eters and possible correlation is discussed and finally, we 
summarize the results in Section 4. 



set of sources has been selected according with the following 
selection criteria: 

• sources classified as blazars, optically identified (e.g. 
iMasetti et ai1l2004l200sh 

• sources also observed in the soft X-ray energy domain 
(0.2 - 10 keV); 

Two further sources 3C 273 and 4C 04.42 have been ex- 
cluded since widely discus sed elsewhere (|Chernvakova et al] 
120071 : iDe Rosa et al.ll200sh . 

Our sample of INTEGRAL selected blazars, is com- 
posed of 9 sources: 7 FSRQ objects with 0.53 < z < 3.67 
and 2 BL Lac objects with z ~ 0.07 and 0.09, but only five 
out of seven FSRQ were available at the time of our analy- 
sis. For the 2 other objects we adopted the hard X-rays BAT 
data. 

Relevant information on the sample are reported in Ta- 
ble 1 where we list the source name, type, optical coordi- 
nates, redshift, Galactic absorption along the line of sight 
according to Dickey & Lockman (1990). In the last column 
of Table 1 we also list the IBIS 20 - 100 keV flux as reported 
by Bird et al. (2007) and the BAT 14 - 195 keV flux as re- 
ported by AGN Catalog of the first 9 months available at 
the time we started this project 0. 

We remark that source with ICR name (IGR 
J22517+2218) is discovered in hard X-ray with INTEGRAL. 
It is the farthest object so far detected by INTEGRAL whose 
nature w as determined a pos teriori through optical spec- 
troscopy (jBassani et al.l l2007). It is worth to note that two 
sources (FSRQs 3C 279 and PKS 1830-211) have been de- 
tected in the Gamma Ray band with the Astro-rivelatore 
Gamma a Immagini L e ggero (AGILE) R (jGiuliani et al] 
120091 ; iTavani et al]|200Sl ; IPittori et all 120091 ) and the Large 



Area Telescope on board the Fermi Gamma-ray Large Area 
Space Telescope (GLAST) 0. The latter satellite also re- 
vealed a third source belo nging to our sample, the object 
BL Lac l|Abdo et. alll2009r i. 



2.2 INTEGRAL DATA 

The INTEGRAL data presente d here are based on point- 
ings with the IBIS instrument ljUbertini et al] 12003). col- 
lected over the period from end of 2002 up to April 2006 
revolution 12 up to 429). Images from the ISGRI detector 
Lebrun et al.ll2003l ) for each pointing have been genereated 
in differe nt energy bands using off-line scientific Analysis 
Software (|Goldwurm et al] 120031 ) OSA version 5.1. Count 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 THE HARD-X RAYS SELECTED SAMPLE 

The INTEGRAL/IBIS total sample is derived from the third 
IBIS/ISGRI survey catalog that included 421 sources for a 
total exposure time of 40Ms; 131 sources out of 421 have 
been identified as AGN. The survey input data set consists 
of all pointing data available at the end of 2006 May, from 
revolutions 12-429 inclusive, covering the time period from 
launch (17 October 2002) to the end of 2006 April. Details 
about this survey are presented by Bird et al. (2007). A sub- 



1 The BAT data have been taken from the on-line archive at: 
http://swift.gsfc.nasa.gov/docs/swift/results/bs9mon For the 
BAT spectra (related to PKS 2149-306 and 0537) we stress that 
in the 22-months survey now available, the 15-45 keV fluxes are 
doubled with respect to the 9-months survey. This effect is taken 
into account with the cross calibration constant that has been 
added when dealing with the fit between soft-X and hard-X rays 
spectra. 

2 An Italian Space Agency (ASI) mission launched on 23 April 
2007 with a key scientific project: the Gamma-Ray observations 
of blazars. 

3 An international and multi-agency mission launched on 11 June 
2008; one of its major scientific goals is to provide new data on 
the Gamma-ray activity of AGN. 
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Table 1. Data for our sample of blazars observed with INTEGRAL, Swift and XMM-Newton. 



Source 


Type 


Broad — Band 




R.A. 


Dec 


redshift 


N Gal 
H 

[10 22 cm~ 2 ] 


j?IBIS 
20-100fc 

, nBAT 

14— 195fc 
[10-"er a .cm- 


eV 
eV> 


1™'-'™"' 
20 — lOOfceV 

rest-frame 
^ 14— 195fceV 

[lO^ers, ■ 


QSO B0836 + 710 


FSRQ 


XMM+INTEGRAL 


08 


41 


24.37 


+ 70 53 42.2 


2.172 


0.030 


6 




63.1 


1ES 0033+595 


BL Lac 


XMM+INTEGRAL 


00 


35 


52.63 


+ 59 50 04.6 


0.086 


0.427 


2 




0.1 


PKS 0537-286 


FSRQ 


XMM+BAT 


05 


39 


54.28 


-28 39 55.9 


3.104 


0.020 


(3) 




(31.6) 


PKS 2149-307 


FSRQ 


XMM+BAT 


21 


51 


55.52 


-30 27 53.7 


2.345 


0.020 


(5) 




(39.8) 


Swift J1656. 3-3302 


FSRQ 


XRT+INTEGRAL 


16 


56 


16.56 


-33 02 09.3 


2.4 


0.220 


2 




50.1 


IGR J22517+2218 


FSRQ 


XRT+INTEGRAL 


22 


51 


53.50 


+ 22 17 37.3 


3.668 


0.050 


4 




100.0 


Bl Lac 


BL Lac 


XRT+INTEGRAL 


22 


02 


42.72 


+ 42 17 16.8 


0.069 


0.022 


3 




0.1 


PKS 1830-211 


FSRQ 


XMM+INTEGRAL 


18 


33 


39.89 


-21 03 39.8 


2.507 


0.260 


5 




100.0 


3C 279 


FSRQ 


XRT+INTEGRAL 


12 


56 


11.17 


+ 05 47 21.5 


0.536 


0.020 


2 
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rates at the source position have been extracted from in- 
dividual images to provide light curves in different energy 
bands. From light curves the average fluxes have been de- 
rived and com bined to produce an average source spectrum 
(for details see lBird et al]|2007r i in the 20 - 100 keV band. 

2.3 XMM-NEWTON DATA 

XMM data of the six blazars (1ES 0033+595, 4C 04.42, PKS 
1830-211, QSO B0836-710, PKS 0537-286 and PKS 2149- 
307) are a combination of proprietary data and public obser- 
vations obtained from the XMM-Newton Science Archivfl 
The raw EPIC Observation Data Files (ODFs) were ob- 
tained from the XMM Science Archive and reduced using 
the standard Science Analysis System (SAS) software pack- 
age (v. 7. 1.0) and the most recent calibration files available 
at the time of the data reduction. We used the EMCHAIN 
and EPCHAIN task for the pipeline processing of the ODFs 
to generate the corresponding event files. The spectra were 
created using X-ray events of pattern 0-12 for MOS and 0-4 
for PN. The source counts were extracted from a circular re- 
gion centred on the source with a radius of 20-40 arcsec and 
the background was derived from two nearby source-free cir- 
cular regions of the same size. Spectra were re-binned using 
GRPPHA to have a minimum of 20 counts in each bin0, so 
that the \ 2 statistic could reliably be used. Details of the 
XMM-Newton observations and the source observed count 
rates are summarized in Table 2. 

2.4 SWIFT/XRT DATA 

Swift /XRT data reduction of four blazars (3C 279, BL Lac, 
IGR J22517+2218, Swift J1656.3-3302) was performed using 
the tool XSELECT v. 2.4. Events for spectral analysis were 
extracted within a circular region of radius 20 arcsec centred 
on the source position. The background was extracted from 

4 http:/ /xmm.esac.esa.int/xsa/index.shtml. 

5 10 counts are the minimum number of counts required to use 
the x 2 minimization technique. 



a circular region with the same radius and located far off 
the source. In all cases, the spectra were binned using GRP- 
PHA. We used Ancillary Response Files (ARFs) and Re- 
sponse Matrix Files (RMFs) available for download from the 
HEASARC Calibration Database (caldb) calibration files at: 
http: / /heasarc. gsfc.nasa.gov /docs /heasarc / caldb / caldb_intro.html 
Table 3 lists the observation ID, the observation start time, 
the integration time and the observed count rate. 



3 BROAD-BAND SPECTRAL ANALYSIS 

The spectral analysis was performed by combining the dif- 
ferent data from each instrument using the X-r ay SPECtral 
fitting software XSPEC v. 11.3.2 (|Arnaudlll996T ). The errors 
are quoted at 90% confidence level for one interesting pa- 
rameter. The energy values are always given in the observer 
frame, if not otherwise specified. Throughout this paper, a 
WMAP (Wilkinson Microwave Anisotropy Probe) Cosmol- 
ogy of H = 70 km s _1 Mpc -1 , Q x = 0.73 and Q m = 1-fi* 
is assumed. 

We have combined XMM spectrum with the IBIS one 
for the blazars 1ES 0033+595, 4C 04.42, PKS 1830-211, 
QSO B0836-710; XRT spectrum with the IBIS one for 3C 
279, BL Lac, IGR J22517+2218, Swift J1656. 3-3302 and 
XMM spectrum with the BAT one for PKS 0537-286, PKS 
2149-307 (see Table 1). 

To account for possible cross-calibration mismatches be- 
tween different instruments as well as to take into account 
flux variations between the observing periods, a multiplica- 
tive constant factor (C) has always been added to the fit. In 
general, the calibration between the various satellites is per- 
formed analyzing the Crab spectr um. In particular, various 
studies of the Crab spectrum (e.g. iKirsch et al"1l2004l , 120051 ) 
have shown that the inter-calibration XMM/INTEGRAL 
and Swift/INTEGRAL is close to 1 (within a few percent). 
We remark that, since XMM, Swift and INTEGRAL obser- 
vations were not simultaneous, flux variations are also pos- 
sible according to the blazar-nature of objects. Being the 
inter-calibration constant C ~ 1, we suggest that, a differ- 
ent value for this constant is likely due to flux variability. In 
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Table 2. XMM-Ncwton observations log and source observed count rates. 



Source 



Obs.ID 



Date of obs. 



Exp. time 
M0S1 

(*) 



Exp. time 
MOS2 

(*) 



Exp. time 
PN 

w 



Counts/s 
MOS1 
(0.2 - WkeV) 



Counts / s 
MOS2 
(0.2 - lOkeV) 
(s- 1 ) 



Counts/s 
PN 

(0.2 - WkeV) 
is' 1 ) 



1ES 0033+595 0094381301 2003-02-01 1204 1535 4209 

PKS 1830-211 0204580301 2004-03-24 30445 30444 20906 

QSO B0836-710 0112620101 2001-04-12 35614 35620 29035 

PKS 0537-286 0206350101 2005-03-20 80610 80765 64661 

PKS 2149-307 0103060401 2001-05-01 23959 23960 20287 



1.39 ± 0.04 
0.58 ± 0.01 
3.88 ±0.01 



1.39 ±0.04 
0.58 ±0.01 
3.85 ±0.01 



3.61 ±0.03 

I. 65 ±0.01 

II. 97 ±0.02 



0.205 ± 0.002 0.226 ± 0.002 0.746 ± 0.004 
0.63 ±0.01 0.62 ±0.01 1.67 ±0.01 



Table 3. Swift/XRT observations log and source observed count rates. 



Source 



Obs.ID Date of obs. Exp. time Counts/s 

XRT XRT 

(0.2 - WkeV) 
(s) (s-i) 



BL Lac 00090042010 2008-08-29 5806 0.20 ± 0.01 

ICR J22517+2218 00036660001 2007-05-26 8300 0.054 ± 0.003 

Swift J1656.3-3302 00035272002 2006-06-13 4785 0.075 ± 0.004 

3C 279 00035019009 2008-11-26 22580 0.198 ±0.003 



our analysis, this constant C between the two instruments 
at low and high energies has been allowed to vary freely 
in the fits. In particular, we have considered the MOS1 and 
XRT instruments as reference detectors in the case of XMM- 
INTEGRAL or BAT data and XRT-INTEGRAL data, 
respectively. Therefore, the inter-calibration constants in 
our fitting procedure are MOS2/MOS1, PN/MOS1, INTE- 
GRAL/MOS1, BAT/MOS1 and INTEGRAL /XRT. How- 
ever MOS2/MOS1 and PN/MOS1 are always ~ 1 within 
few percent. 

3.1 SPECTRAL FITTING RESULTS 

To reproduce the shape of the intrinsic continuum of our 
sources we fitted the combined spectra (0.2 — 100 keV) of 
each source with several models. 

Firstly, we fit the spectra of our sources in the 3 — 100 
keV energy band with the model A, that corresponds to a 
simple power-law absorbed by gas of the Milky Way. The 
results of the fit are summarized in Table 4. The model A 
provides a good fit for all sources. When extrapolating this 



model to the low energy range, a very poor fit in all cases 
but one (PKS 2149-307) was obtained; in fact the data to 
model ratio (see Figure 1) shows a clear systematic residuals 
below ~ 3 keV in the observer frame. From the ratios shown 
in Figure 1, we can see that the spectra of the majority 
of the sources clearly show a deficit of soft X-ray counts 
with respect to the simple power-law model, indicating the 
possible presence of excess absorption and/or of an energy 
break in the intrinsic continuum. 

In order to properly reproduce such a feature we tested 
the impact of different models. 

In the model B we added an extra-absorption in the 
rest-frames of the sources (phabs*zphabs*zpo model in 
XSPEC) and results are presented in Table 5. 

All sources of our sample, with the exception of PKS 
2149-307 and PKS 1830-211, were well fitted with the model 
B, while for the source PKS 2149-307 we have obtained an 
upper limit value on the intrinsic absorbing column density 
and for PKS 1830-211 only a poor quality fit to the data has 
been derived. Details on this source will be discussed later 
in this section. 
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1ES 01133+595 



IGR J22517+2218 



PKS 2149-307 




channel energy (keV) 



BLLac 



channel energy (keV) 



PKS 1830-211 



channel energy (keV) 



Swift J1656-3302 





channel energy (keV) 



channel energy (keV) 



channel energy (keV) 



Figure 1. The data to model A ratio as obtained by fitting a simple power-law over the observed frame 3 — 100 keV energy range, 
and extrapolating over the 0.2 — 100 keV range. The black, red and green bins up to 10 keV are the data of MOS1, MOS2 and PN 
instruments, respectively. 



Then, we tested a broken power-law continuum model C 
(phabs*bknpo model in XSPEC) for the whole sample. Re- 
sults are shown in the Table 6. This model did not improve 
the fit further in most of the cases and, also in those cases 
for which a slightly better \ 2 is obtained, the null hypoth- 
esis probability increases only a few per cent with respect 
to model B. Although strong statistical evidence does not 
exist for a broken power-law respect to the absorbed one, 
however the presence of an intrinsic spectral break cannot 
be ruled out. 



For PKS 1830-211, both models B and C fail to repro- 
duce source data. PKS 1830-211 is a gravitationally lensed 
galaxy and the excess absorption can be attributed to the 
intervening galaxy at redshift z = 0.89. The best-fit of the 
observed spectrum requires both an absorber and a broken 
power-law (model D), confirming the results reported by 
I Zhang et all (|2008h with three XMM-Newton observations 



Obs.ID 0204580201, 0204580301 and 02045804010 and the 
INTEGRAL average data (up to April 29, 2006). In fact, in 
addition to the excess cold absorption, PKS 1830-211 also 
exhibits an energy break at about 4 keV in the observer 
frame. Moreover, our result for the absorption in excess to 
the Gala ctic one, is in a g reeme nt with that obtained previ- 
ously by |Pe Rosa et "all (|2005l l with Chandra observations 
and a subset of the actual IBIS data. Table 7 reports the 
results of our best-fit model. 

Finally, we also checked a Compton reflection compo- 
nent with a powe r-law reflected from neutral ma tter (pexrav 
model in XSPEC. [Magdziarz fc Zdziarskil I1995I D; the reflec- 
tion component is parameterized in terms of R=Q/2n, the 
solid angle in units of 2-7T subtended by the reflecting mat- 
ter, assumed to be observed face-on). This model resulted 
in a worst fit and upper limit values on R for all absorbed 
sources (see Table 8). 



6 In our analysis we have chosen the observation with longer ex- 
posure time. 
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In addition, we have searched for the presence of an iron 
line emission in the spectra of our sources. To this aim, we 
adopted the following method: we exclude the Fe K band, 
fit a continuum, then freeze the continuum, after we include 
the Fe K band and perform the fit again. The comparison 
between the values of the x 2 of the fits performed including 
and excluding the Fe K band provides us a clear indication 
for the real presence of this feature of the observed spectra. 

For all sources, but one (IGR J22517+2218), the appli- 
cation of this method does not produce any clear evidence 
for the presence of this feature in the spectrum. The main 
results for the source IGR J22517+2218 can be summarized 
as follows: the first step described before resulted in a pho- 
ton index r=1.5 and an upper limit on iVfr=5xl0 22 cm -2 
with x 2 /dof=23/29; in the second step we included the iron 
channels and then we added a narrow line (cr=0.01 keV) 
and the fit has been repeated obtaining x 2 /dof=37/44 with 
a probability PFtest > 99.7% implying that this feature is 
real with an Equivalent Width EW=80 ± 50 eV; repeating 
the fit with a free photon index a value of T = 1.6 ± 0.2 
is obtained with an iron line that is significant at the 99%. 
In summary, for this source there is the evidence for the 
presence of a weak line significant at the 99% level. 

As a conclusion, we firmly asses that the an iron emis- 
sion line is not request by our data, with the exception of 
the source IGR J22517+2218. 



3.2 DISCUSSION 

Our broad-band spectral analysis has shown that a deficit 
of photons in the soft X-ray energy domain is a common 
feature in our sample. This phenomenon has bee n previously 
observed in several sampl e s of RL ob jects ( e.g. iPage et al.l 
120051 : iGalbiati et al.ll2005l ; lYuan et alj|2006h and its origin 
is not yet well understood. It can be produced by either 
an intrinsically curved continuum or an extra absorption 
component. 

The broken continuum hypothesis is linked to the 
blazar-type source. In fact, in the hypothesis of an intrinsi- 
cally curved blazar continuum, one can assume that the rela- 
tivistic electrons in the jet follow a broken power-law energy 
distribution with a low energy cut-off, producing a flattening 
in the observed spectrum. This is a purely phenomenolog- 
ical form assumed to reproduce the observed shape of the 
blazar SEDs, without any specific assumption on the ac- 
celeration/cooling me chanism acting on the particles (e.g. 
iTavecchio et al.ll2007I ). Alternatively, it can be also assumed 
there is an inefficient radiative cooling of lower energy elec- 
trons of the jet producing a few S ynchrotron photons to be 
scatt ered at higher energies (e.g. ISikora. Begelman fe Reea 
1 19941 ). 

We emphasize that on basis of our spectral analysis an 
intrinsic spectral break cannot be ruled out (see the val- 
ues of the reduced x 2 and of the null hypothesis probability 
listed in the Tables 5 and 6). In this context, we just noted 
that the existence of a break in the intrinsic continuum is 
tightly related to the electron distribution, but important 
constraints on the curved electron distribution can be de- 
rived only analyzing both Synchrotron and EC components 
in the SED, that is out of the aim of present work. There- 
fore we will focus the following discussion on the absorption 



The absorption explanation takes into account an ab- 
sorbing gas in excess of the Galactic one; in some cases, the 
absorbing medium can be cold (not ionized) and located 
either between the observer and the source along the line 
of si ght at a different redshift with re s pect to the blazar 



or sig nt at a dirterent redsmit witn re s pect to tne blazar 
(e.g.lCappi et al.ll 19971 : iFiore et al.lll998l : IFabian et al"1l200ll : 
iMasetti et al. 20041). or at the same re dshift of the source 
(e.g. IPage et al. I I2OO5I : lYuan et al.l l ioOrj ') and it is namely an 
intrinsic cold absorber. 

For an absorber located at several redshifts between 
the observer and the source, the da mped Ly-alpha sy stems 
are the more plausible cause (e.g. lElvis et al.1 1 19941 ') and 
the line of sight would have to pass through two or more 
very high colum n density damped Ly-alpha systems (e.g. 
iFall fe Pei|[l995l). However, th i s is a rare occurrence (e.g. 
l O 'Flaherty fe Jakobsenl Il997l ; IZwaan. Verheiien fe Briggsl 
1999) and favours the hypothesis of an intrinsic absorber 
with column densities of the order 10 22 — 10 23 cm -2 . An 
intrinsic origin is also confirmed by the soft X-ray spec- 
tral flattening detected so far only in radi o loud objects 
jBechtold et al.l Il994l : iReimers etall 1 19951 : ISiebert et ail 
Il996h . 

The intrinsic hypothesis also suggests that the medium 
could be "warm" at the same stage. In this case, the absorb- 
ing medium, at the same redshift of the AGN (probably, the 
inner part of the dense interstellar medium of the young host 
galaxy), is ionized as due to the pro ximity from the jet (e.g. 
IFabian et alll200ll : lYuan et al.ll2006l ). 

For our sample, the fitting procedure supports the pres- 
ence of an intrinsic (local to the AGN) cold absorber for 
all sources. Nevertheless a warm absorber model (absori in 
XSPEC) can be excluded since the model with an ionized- 
gas produces a worse fit and an inconsistent value for the 
ionization parameter (£ ~ erg cm s _1 and £ > 500 erg cm 
s , values indicating a neutral absorber and a completely 
ionized absorber, respectively). 

The information on the absorption is of crucial impor- 
tance to understand the blazar environment and its interac- 
tion with the jet. 

Absorption originating from the material present in the 
AGN environment seems to be a more convincing explana- 
tion for many sources belonging t o the class of RL AGN 
jPage et al.ll2005l : lYuan et al.ll2006l . and references therein). 
There is some evidence not only of absorption but also of a 
correlation between absorption and redshift, with the more 
distant sources being more absorbed and the increase i n Nh 
with z occurring starting at z ~ 2 (|Yuan et al . 2006). We 
note that more than half of the sources in our sample ap- 
pear in agreement with this scenario, since this fraction of 
objects has a redshift larger than 2. 

As far it concerns the nature of this absorber, we suggest 
it could be associated to the presence of gas in the host 
galaxy illuminated by the blazar jet. 

The excess of intrinsic absorption in the rest-frame 
of the RL objects compared with the lack of absorp- 
tion a bove the Galactic value of their RQ counter- 
parts dCanizares fe White! fl989l; lLawson fe Turner! [l997j ; 
Yuan et alJI 19981: IPage et al.ll200a. 120051 : IVignali et al.ll2003l . 



20051 : iGrupe et al.l "l2006). suggests that an intrinsic differ- 



ence could exist in the local environments of the two classes 
of objects. Since the X-ray spectra of RL AGN appear to 
be dominated by the emission of the jet component, most 
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Table 4. Best fit parameters of a simple power-law model in 3-100 keV energy range. 
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Model A 

Name 1 F o b(3-10keV) 1 F ob(20-l00keV) 3r 4 C 5 Xr/ dof 6 Pnull 



QSO B0836+710 


4.1 


5.2 


1.34 ±0.02 


0.28 ± 0.03 


0.91/1521 


0.996 


1ES 0033+595 


1.1 


1.5 


9 ,- 9 +0.18 


4 2+ 2 ' 5 


0.79/53 


0.865 


PKS 0537-286 


0.5 


2.3 


1.19 + 0.06 


0.6 +.0.2 


1.05/329 


0.253 


PKS 2149-307 


0.6 


4.9 


1.41 + 0.06 


1.9 ±0.5 


0.91/225 


0.835 


PKS 1830-211 


1.4 


4.4 


1.22 + 0.04 


0.5 ±0.1 


1.03/442 


0.315 


Swift J1656.3-3302 


0.5 


2.1 


1 60+ ' 33 


1 6+ 20 


1.02/18 


0.434 


IGR J22517+2218 


0.3 


3.8 


i 4 o+0.59 




0.66/14 


0.817 


Bl Lac 


0.6 


2.3 


2 01+ ' 43 
Z - Ui -0.41 


4.3+^ 


0.59/11 


0.840 


3C 279 


0.5 


2.1 


2 00+ ' 24 
Z.UU_ 23 




1.06/19 


0.382 



1 Observed flux in the 2 — 10 keV (or 20 — 100 keV) energy range in units of 10 11 erg cm 2 s 1 . 3 Photon index, related to the 
spectral index a (where F v oc u^ a ) by a = T - 1. 4 IBIS/MOS, BAT/MOS or IBIS/XRT cross-calibration constant. 5 Reduced 
chi-squared to degrees of freedom. 6 Null hypothesis probability that is the probability of getting a value of x 2 as large or larger than 
observed if the model is correct. If this probability is small then the model is not a good fit to the data. 



Table 5. Best fit parameters of an absorbed power-law model in 0.2-100 keV energy range. 
Model B 

Name 1 F b(2-wkeV) 1 Kb(2o-iookeV) 2 N H 3 r 4 C 5 Xr/dof 6 P nu u 



QSO B0836+710 


4.0 


4.9 


0.07 ±0.02 


1.353 ±0.005 


0.28 ± 0.03 


1.02/2421 


0.206 


1ES 0033+595 


1.2 


1.5 


0.20 ±0.02 


2.47 ±0.05 


3.5 ±0.8 


0.81/229 


0.985 


PKS 0537-286 


0.5 


2.3 


0.30 ±0.05 


1.22 ±0.01 


0.7 ±0.2 


1.03/972 


0.235 


PKS 2149-307 


0.6 


4.8 


< 0.05 


1.45 ±0.01 


2.1 ±0.5 


0.99/849 


0.543 


Swift J1656. 3-3302 


0.5 


2.0 


7 o+s.i 
' - a -3.8 


1 66+ ' 22 


1 9+ 1 ' 7 


1.03/42 


0.422 


IGR J22517+2218 


0.2 


3.5 


3.6±H 


, KK +0.21 
L66 -0.19 


7 1+6.7 
'■ L -3A 


0.85/49 


0.770 


Bl Lac 


0.5 


2.3 


0.20 + 0.07 


2 01+ ' 17 




0.81/31 


0.762 


3C 279 


0.5 


2.2 


0.04 + 0.02 


1.75 ±0.06 


2 5+°' 7 
Z -°-0.6 


1.07/79 


0.311 



1 Observed flux in the 2 — 10 keV (or 20 — 100 keV) energy range in units of 10 -11 erg cm~ 2 s _1 . 2 Intrinsic column density of 
hydrogen in units of 10 22 cm~ 2 . 3 Photon index, related to the spectral index a (where F v oc v~ a ) by a = T — 1. 4 IBIS/MOS, 
BAT/MOS or IBIS/XRT cross-calibration constant. 5 Reduced chi-squared to degrees of freedom. 6 Null hypothesis probability that is 
the probability of getting a value of \ 2 as large or larger than observed if the model is correct. If this probability is small then the 

model is not a good fit to the data. 
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Table 6. Best fit parameters of a broken power-law model in 0.2-100 keV energy range. 



Model C 
Name 



1 F o b{2-l0keV) 1 F ob(2()-10akeV) 7 ri,r2 



4 C 5 x 2 r /doi 6 P null 



QSO B0836+710 



1ES 0033+595 



PKS 0537-286 



PKS 2149-307 



IGR J22517+2218 



Bl Lac 



4.0 



1.2 



0.5 



0.6 



Swift J1656. 3-3302 0.6 



3C 279 



0.2 



0.5 



0.5 



4.9 



1.6 



2.3 



4.8 



2.0 



3.5 



2.2 



2.2 



1 26+ 01 
1.354 ± 0.005 

, qo+0.16 
t-38_ 13 

r, o 7 + 0.06 

*- a '-0.05 

0. 53t° f 2 
1.21 ± 0.01 

1 47+0.01 
1 - f " -0.02 
1 41 +0.04 

] n3+ ' 17 

1 96 +0 ' 30 

1. »o_ 27 

rj co+0.49 
u - oo _0.69 
1 69+ ' 23 

1 06+ ' 57 

-0.73 

2 02+ ' 34 
z - uz -0.12 

1 46+°' 10 
1.78 ±0.06 



0.98 ±0.08 0.28 ±0.03 



1 39+ 012 



0.72 ± 0.09 



9 7+4.5 
z -'-l.l 



o ,+45.5 
8 - d -3.4 



1 19+0-27 
1 - ± "-0.26 



i 71 +0.75 
1 -' 1 -0.1S 



1 16+ 027 
± - 1D -0.20 



9 S+0' 7 



0.7 ±0.2 
1.9 ±0.5 



1 S + 2 ' 2 
1.8_ 1 4 



7.8«; 7 



4 6+1-7 
9 7+0.8 



1.02/2420 0.274 

0.85/228 0.947 

1.03/971 0.282 

0.99/848 0.579 

0.87/41 0.708 

0.82/48 0.803 

0.89/30 0.633 

1.01/78 0.446 



1 Observed flux in the 2-10 keV (or 20-100 keV) energy range in units of lO-^crgcm^s- 1 . 4 IBIS/MOS, BAT/MOS or IBIS/XRT 
cross-calibration constant. 5 Reduced chi-squared to degrees of freedom. 6 Null hypothesis probability that is the probability of getting 
a value of x 2 as large or larger than observed if the model is correct. If this probability is small then the model is not a good fit to the 
data. 7 Spectral index below the break (fi) and above the break (f^). 8 Observed break energy in units of keV. 



Table 7. Best fit parameters of a broken power-law model with 


an absorption 


in excess in 


0.2-100 keV 


energy range. 


Model D 










Nam e 1 Fob(2-10keV) 1 F ob(20- WOkeV) 2 ^H 


7 ri,r 2 




4 C 


5 X 2 /dof 6 P null 


PKS 1830-211 1.4 4.3 1 - 5 -oi 


0.95 ±0.05 

1 qo + 0.10 


4-0ll^ 2 4 


8 +0 - 3 

u -0.2 


1.06/856 0.101 



1 Observed flux in the 2 — 10 keV (or 20 — 100 keV) energy range in units of 10 -11 erg cm~ 2 s — 1 . 2 Intrinsic column density of 
hydrogen in units of 10 22 cm — 2 . 4 IBIS/MOS, BAT/MOS or IBIS/XRT cross-calibration constant. 5 Reduced chi-squared to degrees of 
freedom. 6 Null hypothesis probability that is the probability of getting a value of x 2 as large or larger than observed if the model is 
correct. If this probability is small then the model is not a good fit to the data. 7 Spectral index below the break (ri) and above the 

break (r2). 8 Observed break energy in units of keV. 



Table 8. Upper limits to R under the assumption of a pexrav model. 



Name 


QSO B0836+710 


1ES 0033 + 595 


PKS 0537-286 


PKS 2149-307 


PKS 1830-211 


Swift J1656. 3-3302 


IGR J22517+2218 


Bl Lac 


3C 279 


it 








0.09 


0.08 





0.5 


3E-6 





3E-4 
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likely the absorbing material is linked to the presence of this 
component. This working scenar io is supported by a theory 
presented, e.g., in iFerraril (|l998h . in which the confinement 
of the radio jet needs the presence of a cold medium (local 
to the AGN) so that jet expansion is suppressed. 

In our sample we have found neither an evidence 
of an iron emission line at 6.4 keV nor evidence 
of a Compton reflection "hump" above 10 keV. The 
lack of these spectral features - interpreted as the re- 
sult of reprocessing of the primary continuum by cold 
matter aroun d the X-ray source, pre s umably the ac- 
cretion disk (|Lightman fc Wlift d Il988l: iGuilbert fc ReesI 
19881; iGeorge fc Fabianlll99ll; ISambruna fc Eracleous! 1199a 
Sambruna. Eracleous fc Mushotzkvl I1999T ) - can be ex- 



N distribution 



plained either with the out-and-out absence of reprocess- 
ing features in blazar-type objects or with the hypothesis 
that these features are difficult to be detect in strong radio 



sources being likely dilute by the jet \ Sambruna ct al. 2006). 

There is also evidence that the spectral shape does not 
change between 0.2—10 keV and 20 — 100 keV energy ranges, 
supporting the idea that, above 10 keV, no additional com- 
ponent to the power-law (e.g. Compton reflection hump) is 
present in the spectrum. 

The value of the cross-calibration constant (see Tables 
B and C) shows that a strong flux variability is present 
in QSO B0836+710, 1ES 0033+595, PKS 1830-211, IGR 
J22517+2218, BL Lac, and 3C 279. In particular, we derived 
that IGR J22517+2218 showed severe intensity changes be- 
tween the INTEGRAL and the Swift observations with a 
value of the XRT/IBIS inter-calibration factor of about 7. 
This is not a surprise because of the nature of the object 
and data being not simultaneous. In addition, we remark 
that INTEGRAL/IBIS data set correspond to an average 
performed on observations spanning on 3 years, while the 
data for the low-energy range (XMM and XRT) correspond 
to observations with a time baseline of few hours. 



3.3 PARAMETER DISTRIBUTIONS AND 
CORRELATIONS 

Although the number of sources of our sample is not large, 
we have also analyzed the distribution of the spectral pa- 
rameters (absorbing column density, spectral index, red- 
shift) and their possible correlations. This type of analy- 
sis has been perfo rmed, so far, only between 0.2 — 10 keV 
l)Yuan et al.ll2006h . Our study is still preliminary and has 
been performed with a detailed comparison among the re- 
sults of our own spectral analysis and those already available 
in literature in a limited energy range. The reference sample 
includes 3 5 RL objects - with 19/35 sources i dentified as 
blaza rs - (|Galbiati et al.ll2005l : IPage et al.ll2005l : lYuan et all 
l2006h and shows an absorption of unknown intrinsic nature. 
In order to avoid any problem of "orientation dependence 
bias" when studying the Nh and T distribution, we have se- 
lected, from the quoted sample, only the blazar-type objects. 
In particular, for the Nh study we have accounted for the 
blazars with a specific value of the intrinsic absorption (we 
have not accounted for the upper limits on Nh), obtaining 
a sub-sample of 16 objects. 

In the framework of a possible A r ij-redshift relation, 
we also combined our results with the larger s a mple of RL 
objects (jGalbiati et all 120051 ; IPage et al.l 120051 : lYuan et all 




Figure 2. Distribution of the X-ray absorption column density 
(Nh) for our sample objects over the 0.2 — 100 keV energy range. 
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Figure 3. The Nh distribution of a larger sample of 16 blazars 
over the 0.2 — 10 keV energy range. See text for more details. 



2006), counting out three sources belonging to our sam- 
ple (QSO B0836+710 PKS 2149-307 and PKS 0537-286). In 
fact, it is interesting to check whether the 9 selected blazars 
in the sample and analyzed in a more wide energy band, fol- 
low the trend of N h with the redshift previously found by 
lYuan et alj l|2006l ). If this would be the case, we would have 
further evidence that the absorbing column density, in these 
sources, is linked to the distance and thus a confirmation of 
a cosmic evolution effect. 



3. 3. 1 ABSORPTION DISTRIB UTION 

We made a detailed analysis of the rest-frame hydrogen col- 
umn density for our selected blazars, in conjunction with 
the intervening absorption in excess to the galactic one. 

The Nh distribution shown in Figure 2 appears to be 
asymmetric, with a peak located at about log(Nn) = 21 
and tails extending towards low and high values up to 
log(Nn) = 20 and log(Nn) = 23, respectively. It is worth to 
note that the statistical significance of this result is strongly 
affected by the limited number of objects in our sample. 
Therefore we compare our sample with the larger sub-sample 
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Figure 4. The X-ray absorption column density Njj as a function 
of the redshift for our sample. 
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Figure 5. The X-ray absorption column density Nh as a function 
of the redshift, combining the results of our sources (filled rhombs) 
and the ones of a larger sample (plus and cross signs) analyzed 
in the soft X-ray energy range. 



from the literature observed in the 0.2 — 10 keV energy range 
from XMM-Newton that includes 16 RL objects identified 
as blazars and with a specific value of the absorption. We 
underline that this one includes only one object belonging 
to our sample (the blazar QSO B0836+710). The N H dis- 
tribution of the larger sample is shown in Figure 3, with a 
peak located at log(Nn) = 22 and a tail extending towards 
lower values up to log(Nn) = 20. Comparing the distribu- 
tions in Figure 3 and in Figure 2, we note that sources with 
absorbing column density around log(Nn) = 23 are "miss- 
ing" . This evidence is expected in view of the fact that in the 
soft X-ray energy range (0.2 — 10 keV) the sources with an 
absorbing column density Nh ^ 1-3 10 23 cm -2 are almost 
completely obscured. 

This occurrence strongly suggests that data above 10 



keV represent a crucial tool to construct an unbiased selec- 
tion of heavily absorbed objects. 

The Nh -redshift trend of our sample is shown in Figure 

4. 

The co mparison of t his N h-z plot with those previously 
analyzed by |Yuan~et? al. ( 2006) at various redshifts for a sam- 
ple of 32 RL sources, seems to indicate that the property of 
the X-ray absorbing gas evolves with cosmic time. This trend 
is shown in Figure 5. where we added 3 objects analyzed by 
iGalbiati et alj (|2005h . In fact, at first glance, a correlation 
between Nh and redshift seems to exist and, compared to 
low redshifts (z < 2), a shift toward higher values of Nh 
seems to be indicated at high redshifts. 

The sample of 32 sources includes 3 objects of our sam- 
ple (QSO B0836+710, PKS 0537-286 and PKS 2149-307) 
and the Nh values for these 3 sources as obtained through 
our spectral analysis are in agreement with those presented 
in the work of Yuan and collaborators (2006) . 

The existence of a correlation is also confirmed by the 
statistical analysis: we obtain a correlation coefficient (R) 
around 0.7 (Figure 5), when the results correspo nding to our 
sampl e are added to the sou rces i nvestigated by Page et all 
i|2005h . IGalbiati et all (|2005l ) and I Yuan et HI l|2006l ). 

All together these results seem to support the hypoth- 
esis of an intrinsic absorber. Three further hints support its 
intrinsic nature: 

• the existence of a spatial isotropy of the inter-galactic 
absorbers is most unlikely; 

• the X-r ay absorption is associated with RL but not with 
RQ objects ( Yuan fc Brinkmannl fl999). This evidence sup- 
ports the hypothesis that the absorber are physically asso- 
ciated with the RL quasars; 

• the change of the Nh distribution that seems to occur 
at z ~ 2 and the decrease of the slope in the Nh-z relation 
at z > 2 could indicate an effect of cosmic evolution. This 
effect could be due, e.g., to a rate of stellar formation. In 
fact, an increase of the stellar formation at z ^ 2 could 
reduce the amount of gas in the host galaxy. 



3. 3. 2 PRO TON-INDEX DISTRIB UTION 

The photon-index distribution was then investigated by 
adding the FSRQ 3C 273 and 4C 04.42 in our analysis 
(r ~ 1.6 and ~ 1.3, respectively, as obtained by our broad- 
band analysis). The F distribution for the complete sample 
is shown in Figure 6 and the mean value obtained in 0.2— 100 
keV energy band is F = 1.39 ± 0.01. 

For comparison, in Figure 7 we also show t h e pho - 
ton index distribut ion derived from Page et all |2005), 
IGalbiati et al. | i|2005h and I Yuan et"aT1 l|2006h for 35 sources. 
The F distribution of our sample shows a peak at F ~ 1.4, 
whereas the V distribution of the RL objects shows a peak 
corresponding to a steeper value of F ~ 2. However, we note 
here that our sample is mainly composed by FSRQ, while 
the larger sample contains different types of RL objects - 
not all identified as blazars. Therefore, in order to perform 
a meaningful comparison, we have studied the photon index 
distribution for sources - in the larger sample - identified 
as blazars (19 objects). The peak of the F distribution (see 
Figure 8) is, again, around 2. 

To account for the fact that most of the sources in our 
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Figure 6. Distribution of the power-law photon-index for our 
sample in the 0.2 — 100 keV energy range. 



Figure 9. Photon-index distribution for 15 sources identified as 
FSRQ. 



r distribution 




Figure 7. Distribution of the photon-index resul ti ng by the anal- 
ysis p erfor med in the works o f |Paee et alj ||2005| ), iGalbiati et al.l 
[1200511 and lYuan et alj <2006f> and including 35 RL objects. 



r distribution 




Figure 8. Photon-index distribution for 19 sources identified as 
blazars. 



sample are FSRQ-type objects and that FSRQs in X-rays 
show a lower value of V with respect to BL Lac objects, we 
studied the F distribution including only the FSRQ sources 
(15 objects) in the larger sample. In this case, the F distri- 
bution in Figure 9 shows a peak at F ~ 1.6, a value steeper 
than that for our blazars sample. 

The differences between the values of F of our sample 
and those of the larger sample could be due to the hard X-ray 
selection of our INTEGRAL/IBIS sample, which is clearly 
biased towards the natter values of the photon index. 



4 SUMMARY AND CONCLUSIONS 

In this work we presented a broad-band X-ray spectral study 
of a sample of 9 blazars (with a redshift range of 0.1 < z < 
3.7) observed with INTEGRAL, XMM-Newton and Swift. 
The main results can be summarized as follows: 

• The broad-band spectra of all selected sources are well 
reproduced with a power-law model absorbed by an amount 
of gas in excess to the Galactic one (Nh ~ 10 20 -10 23 cm -2 
in rest-frame of the source; only an upper limit of Nh ~ 5 
10 20 cm- 2 has been derived for the FSRQ PKS 2149-307). 

• The absorption seems to be a signature of a cold in- 
trinsic absorber, confirming and extendi ng to larger sample 
previous results quoted in the liter ature (|Cappi et al J 1 19971 : 
iPaee et al.ll2005l : lYuan et alj|2006h . 

• The present work provides a further confirmation of the 
existence of a Nh -redshift trend, obtained for a large sample 
of RL objects (not only blazars). 

• The broad-band analysis of our sample of blazars re- 
vealed a harder spectrum with a photon-index of the order 
of r ~ 1.4, compared to the value obtained with the distri- 
bution including the FSRQ sources of the larger sample (15 
objects). Such a difference could be due to the hard X-ray 
selection of our INTEGRAL/IBIS sample which is clearly 
biased towards flatter values of the photon index. 

• We have found no evidence of reflection components (re- 
flection "hump" and iron emission line), with the exception 
of the source IGR J22517+2218 that shows the presence of a 
weak iron line. This result is expected in blazar-type objects. 
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In conclusion, the analysis presented here has shown 
that our INTEGRAL-sample selection favours objects heav- 
ily absorbed and with a natter value of spectral index. On 
the other hand, the present broad-band analysis of INTE- 
GRAL/IBIS, XMM-Newton and Swift observations of RL 
QSOs confirms that the observed flattening is common in 
these objects, and is clearly detected in eight quasars of 
our sample (8/9). The assumption of an intrinsic origin and 
a cold nature for the absorber is c onsistent w it h pre vious 
results up t o 10 keV obtained by iPage et all (120051 ) and 
lYuan et af] (|200fj ) with XMM-Newton data only. However, 
a broken power-law model, as an alternative explanation for 
the deficit of soft photons observed in the majority of our 
sources, cannot be ruled out by the data. 
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